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Charmonium and bottomonium production is studied using /i + /x — data samples collected by the CDF and DO 
experiments during the 1992-96 pp collider run at the Fermilab Tevatron. The inclusive cross sections as a function 
of the transverse momentum of reconstructed quarkonium states are measured. The results are compared with 
theoretical predictions, which take into account different quarkonium production mechanisms. 



1 Introduction 

During the 1992-96 collider Run I, the Tevatron 
delivered over 150 pb _1 of integrated luminosity. 
The charmonium and bottomonium production 
rate results presented in this paper were measured 
by the CDF and DO collaborations and are based 
on dimuon subsets of the data. In the quarko- 
nium analyses, differential cross sections for J/tp, 
tp(2S) and three T states have been measured us- 
ing their \x + \x~ decay channels. Two muons in the 
final state were used as a trigger signature. 

In early studies of quarkonium production the 
dominant contributions were assumed to come 
from the lowest order Feynman diagrams involv- 
ing gluon fusion — either directly into quarkonium 
states and recoiling partons or through a bb pair 
followed by decays B — > — or, in the case 
of the J/tp or the T, radiative decays of \ c or Xb 
mesons.]] Various methods are used to disentan- 
gle the three sources, providing information about 
charm and bottom production and fragmentation 
mechanisms at low transverse momenta. 

Previous CDF measurements of J /tp and 
tp(2S) production rates during the 1988-89 collider 
rurn showed production cross sections consider- 
ably larger than contemporary theory predicted. 
This drew theoretical interest, but at that time 
the question of whether or not the excess could be 
attributed to a large prompt component was not 
addressed. It has been pointed out that in addi- 
tion to gluon fusion, the gluon fragmentation pro- 
cesses are also important sources of quarkonium 
productionu We briefly summarize the proposed 
improvements in the theoretical description of the 
observed production rates. 



2 Charmonium total cross sections 

Muons identified in the central pseudorapidity 
range, < 0.6, were used in these analyses. Ad- 
ditionally, the DO detector provides muon coverage 
in the forward pseudorapidity region, allowing pro- 
duction rate studies in the range 2.5 < \rj\ < 3.7. 

In the CDF analyses, both muons of the re- 
constructed \& candidate are required to have px 
greater than 2.0 GeV/c, and at least one muon of 
the pair is required to have pt > 2.8 GeV/c. The 
reconstructed J/tp or tp{2S) candidates must have 
p T > 5.0 GeV/c. About 22,000 J/tp events and 
800 tp(2S) events are reconstructed in data sam- 
ples of 15.4 pb _1 and 17.8 pb _1 , respectively. In 
the DO analysis, muons with transverse momen- 
tum greater than 3 GeV/care reconstructed in the 
central and forward pseudorapidity regions, and 
dimuon candidates are required to have px > 8 
GeV/c. About 4,000 $ events in the central and 
about 500 events in the forward rapidity ranges 
are found from a fit to a Gaussian function and 
physics-motivated background in the samples of 
60 pb _1 and 9.3 pb _1 , respectively. The DO track- 
ing momentum resolution is such that it does not 
allow separation of the tp(2S) and J/tp states. 

CDF measures the product of dimuon 
branching ratio times integrated cross sec- 
tion to be 17.35 ± 0.14 ± 2.79 nb for J/tp and 
0.57 ± 0.04 ± 0.09 nb for tp(2S) in the central ra- 
pidity region with px > 5.0 GeV/c, and DO mea- 
sures 1.96 ± 0.16 ± 0.63 nb for J/tp in the central 
region and 0.40 ± 0.04 ± 0.04 nb in the forward re- 
gion with px > 8.0 GeV/c. 

3 Charmonium from b Decay 



"Representing the CDF and DO Collaborations. 
b ^ represents both J/tp and ip(2S) mesons. 



The CDF collaboration, using a silicon vertex de- 
tector, separated the J/tp and tp(2S) samples into 
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their components arising from b decay and from 
prompt production by analyzing the proper decay- 
length (ct) distributions. The ct distribution is 
fitted to three components: an exponential convo- 
luted with a Gaussian resolution function for the 
b hadron decay component, a Gaussian function 
centered at for prompt production, and a Gaus- 
sian function with positive and negative exponen- 
tial tails to describe the background, both combi- 
natorial as well as from sequential b — > pT c — > p + s 
decays. The samples are subdivided into ranges 
of pt(^ + H~) and fitted separately for each range. 
Figure [j] shows the J/ip and ip(2S) fractions from 
b decay as a function of pr- These fractions are 
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Figure 1: Fraction of J/ip and ip(2S) from b decay as a 
function of pp. 
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Figure 2: 6 quark cross section determined from the J/ip 
(top) and ij(2S) (bottom) samples. 



then convoluted with the charmonium px spectra 
to give the b cross sections, as shown in figure |^. 
The results are within a factor of 2-3 of the NLO 
QCD prcdictionp using a central value of the input 
parameters|| and are consistent with other Teva- 
tron b cross section resultsQ 

4 Prompt Charmonium Production 

The cross section for prompt charmonium produc- 
tion is in disagreement with theoretical predic- 
tions based on color-singlet production of bound 

c The input parameters have the following values: renor- 
malization and fragmentation scales /ip = fiji = /iq = 
v/(Wj +Py), the mass of the b quark, mj,=4.75 GeV/c, 
the Peterson fragmentation parameter, e;, = 0.006. 



cc statescl The rate of prompt ip(2S) production is 
about a factor of 50 larger than predictions based 
on such a model. All prompt ijj(2S) are believed 
to be directly produced since \c states with suffi- 
cient mass to decay to ip{2S) lie above the thresh- 
old for strong decays to DD meson pairs. How- 
ever, prompt J/ip are produced not only directly, 
but also via \ c radiative decays. CDF has deter- 
mined the fraction of the prompt J /ip sample com- 
ing from the \c decay by fully reconstructing the 
decay \c — ► J/^Pl- Photon candidates detected in 
the central electromagnetic calorimeter with en- 
ergy greater than 1 GeV and having no charged 
track pointing to the same calorimeter tower are 
combined with p + p~ pairs consistent with the 
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3/ip mass. A peak containing 1230 ± 72 Xc can- 
didates is observed in the mass difference distri- 
bution AI(n + fi~j) — M(n + fj,~). The background 
under the peak has been modelled by embedding 
simulated ir° and 77 decay photons in real J/ip 
events. The fraction of J/ip coming from prompt 
Xc decay, measured for 4 different px bins, ranges 
from about 32% in the 4-6 GeV/c bin to 28% in 
the bin px > 10 GeV/c. The DO measurement of 
this fraction is consistent with these resultsQ Mul- 
tiplying the total prompt J/ip cross section by the 
Xc fraction shows that the rate of J/ip production 
from Xc is within a factor of 2-3 of the theoretical 
prediction, but, as with the ip(2S), the remaining 
direct J/ip cross section is about a factor of 50 
larger than the color-singlet prediction. 

CDF Preliminary 



CDF Preliminary 



> 
O 



L_ 





M( 3 S l , B >) = 1 1 


± 0.8 x 10" 3 GeV 3 


Forcing J/fk 


M('S^, 3 Pf ) ~- 




- and ip' \\ 


= 29.9 x 10' 3 GeV 3 \ 


amplitudes \A 






to have the \ \ 






same ratio \ 




A ftp — > mV" 


— Sum 






: — *Sf , 




\ oN ; 


igfs) 3p(B) \ 






: — Singlet 






- ---- V' Feeddow 
, i.i.i 


, i \ i , 


i , i , r--, i , 



2 4 6 8 10 12 14 16 18 20 

P T (J/V) (GeV/c) 

Figure 3: Prompt direct J/ip sample, compared to the the- 
oretical prediction, with color-octet components fitted si- 
multaneously to the J/ip and ip(2S) distributions. 

One proposal to explain the observed prompt 
charmonium production rates is tp include cc pairs 
produced in a color-octet stateQ The initial pro- 
duction can be calculated perturbatively and can 
be used to predict the pr dependence of the cross 
section. The transition to a color-singlet state 
needed to form a bound cc particle proceeds via 
soft gluon emission. This latter process cannot be 
calculated perturbatively, so the normalization is 
found by fitting the theory to the data. Figures 
H and U show the prompt J/ip and ip(2S) cross 
sections and the corresponding theoretical predic- 
tions when the fitted color-octet contributions are 
includedB The predictions of the model can be 
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Figure 4: Prompt direct tp(2S) sample, compared to the 
theoretical prediction, with color-octet components fitted 
simultaneously to the J/ip and ip(2S) distributions. 




Figure 5: The cross sections - Br cicr /dpx for 2.5 < \rj\ < 
3.7 (left), and Br-da/dr] for p T > 8 GcV/c(right). The solid 
curve represents the expected contribution from b quark 
fragmentation. 



further tested by fixed target hadro- and photo- 
production experiments and by measurement of 
the ip(2S) polarization in pp collisions. 

5 Forward Charmonium Production 

The DO collaboration determined the differential 
cross section, da/dpt, shown in Figure ||, using 
candidates reconstructed in the forward region 
This is the first measurement of the "J 7 cross sec- 
tion at large pseudorapidity. Figure |B| combines 
data from the central and forward $ analyses for 
transverse momentum greater than 8 GeV / c. The 
measurements in the central region are consistent 
with having no r\ dependence. However in the for- 
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Figure 6: X(XS) differential cross section, compared to the 
theoretical prediction. The dotted line shows the color sin- 
glet contribution, while the dashed lines are the color octet 
components fitted to the data. The solid line is the sum of 
all contributions. 

ward region the measured cross section is approx- 
imately a factor of 5 lower than the central region 
values. The data points are compared in Figure 
|5| with a preliminary theoretical prediction of the 
contribution from the b quark decay to 'J'. 

6 Bottomonium Production 

The CDF collaboration has also published^ pro- 
duction cross sections for T(15), T(2S), and 
T(35) states based on a data sample of 16 pb _1 . 
All three states combined yield a total of about 
1800 candidates reconstructed in the /J. + (J.~ de- 
cay channel. The DO collaboration has measured 
combined T cross section based on_events recon- 
structed in a 6.6 pb _1 data samplcllil 

Recently, a theoretical prediction including 
color-octet contributions was fitted teethe T(15) 
and T(2S) differential distributions Figure g 
shows T(15) fit results, which describe the shape 
of the pt distribution well. 

Increased statistics using CDF's full 110 pb -1 
data sample will allow finer binning in pr, improv- 
ing the experimental description of the shape. Re- 
construction of the Xb states via T7 decay — while 
difficult due to the small mass difference between 
Xb and T — would provide an additional probe of 
the underlying production mechanisms. 



7 Conclusion 

Measurements of the differential production cross 
sections for J/tp, ip(2S), and three T states 
have been made at the Tevatron Collider. The 
prompt and 6-decay components of both charmo- 
nium states have been extracted. The prompt J/ip 
cross section has been further subdivided into its 
direct and Xc components. 

These measurements provided the impetus for 
new theoretical models, such as the color-octet 
model, which show potential to explain charmo- 
nium production in pp collisions. Additional ex- 
perimental results, such as measurement of the 
tp(2S) polarization and reconstruction of Xb states 
may provide additional insight into the underlying 
production mechanisms. 
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